Orbital magnetic field effects in spin liquid with spinon Fermi sea: 
Possible application to k-(ET) 2 Cu 2 (CN) 3 
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We consider orbital magnetic field effects in a spin liquid phase of a half-filled triangular lattice 
Hubbard system close to the Mott transition, continuing an earlier exploration of a state with spinon 
Fermi surface. Starting from the Hubbard model and focusing on the insulator side, we derive an 
effective spin Hamiltonian up to four-spin exchanges in the presence of magnetic field, and find that 
the magnetic field couples linearly to spin chirality on the triangles. The latter corresponds to a 
flux of an internal gauge field in a gauge theory description of the spin liquid, and therefore a static 
such internal flux is induced. A quantitative estimate of the effect is obtained using a spinon mean 
field analysis, where we find that this orbital field seen by the spinons is comparable to or even 
larger than the applied field. We further argue that because the stiffness of the emergent internal 
gauge field is very small, such a spinon-gauge system is strongly susceptible at low temperatures to 
an instability of the homogeneous state due to strong Landau level quantization for spinons. This 
instability is reminiscent of the so-called strong magnetic interaction regime in metals with the usual 
electromagnetic field, but we estimate that the corresponding temperature-magnetic field range is 
significantly broader in the spinon-gauge system. 



I. INTRODUCTION 

Experimental studies^Si 3 . of the quasi-two-dimensional 
organic material k-(ET) 2 Cu2(CN) 3 strongly suggest a 
spin liquid state in the insulating phase at ambient pres- 
sure. k-(ET)2Cu2(CN)3 is a strongly correlated half- 
filled Hubbard system on an almost isotropic triangular 
lattice. The material is just near the boundary of the 
insulator-metal transitioni2iiMi£i£i£ It is an insulator at 
ambient pressure with a charge gap of order 350 K, and 
shows no signs of magnetic ordering down to 32 mK de- 
spite a relatively large exchange coupling J ~ 250 K. 
Experiments find that this spin liquid maintains a finite 
susceptibility and a finite 1/(T]T) for the 13 C nuclear 
spin relaxation rate down to low temperatures^^ and 
also that this insulator has large spin entropy^ 3 . 

These observations and related numerical studies^ led 
the present author^ and also the authors of Ref. llOl to 
propose a spin liquid state with spinon Fermi surface as 
a likely candidate. In Ref. 0, this state was shown to 
be energetically favored in a spin model with Heisenberg 
and ring exchanges appropriate for the description of the 
insulator, while Ref. [T(J used the mean field slave parti- 
cle approach and also derived an effective gauge theory 
description of the proposed spin liquid. 

Here, we examine the response of this spin liquid to 
strong magnetic fields. One motivation is to look for pos- 
sible direct probes of the spinon Fermi surface similar to 
the ones used in metals such as magneto-oscillations or 
magnetoacoustic resonance. Another motivation is the 
observations in Refs. H andU of significant broadening 
of the 13 C NMR line in a large magnetic field up to 
9 T at low temperatures below 10 — 30 K. Systematic 
studies 3 reveal that the lines broaden symmetrically and 
the broadening increases with increasing field and also 
with decreasing temperature. At the present, it is not 
clear whether this anomalous behavior is intrinsic, or is 



the result of extrinsic disorder effects in the critical spin 
liquid. 

A simple intuition about the Mott insulator would be 
that the charge motion is suppressed and only spin de- 
grees of freedom couple to the magnetic field. If this were 
the case, then one would not expect any strong intrinsic 
effects in the presence of the magnetic field. In partic- 
ular, in the spin liquid state with spinon Fermi surface 
and just the Zeeman spin coupling to the applied field 
one would expect only Pauli spin paramagnetism. We 
argue, however, that orbital effects need to be carefully 
included when analyzing the response of the spin liquid 
in k-(ET) 2 Cu2(CN) 3 . This is because charge fluctua- 
tions, which become more prominent in the vicinity of 
the insulator-metal transition, also induce an effective 
"orbital" coupling to the magnetic field. The effect of 
such coupling on the proposed spinon Fermi sea state 
is further amplified by the fact that this phase itself is 
stabilized against other competing spin liquids or the an- 
tiferromagnctically ordered state by charge fluctuations 
that produce the four-spin ring exchange termsAiS 

First, we show that in the presence of the external 
magnetic field the spinons effectively experience an "in- 
ternal" orbital field that is comparable to and maybe even 
larger than the applied field. This is despite the fact that 
spinons do not transport electrical charge, but follows 
when we derive an effective spin Hamiltonian from elec- 
tronic degrees of freedom in the presence of the applied 
magnetic field. We find that the magnetic field couples 
linearly to the spin chirality on the elementary triangles. 
The effective description of the spin liquid state contains 
spinons coupled to a dynamically generated "internal" 
gauge field. The physical meaning of the flux of the in- 
ternal gauge field is precisely the spin chirality, and the 
external magnetic field therefore induces a static internal 
flux seen by the spinons that is comparable to the applied 
field. 
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Second, we argue that because the stiffness of the inter- 
nal gauge field is very small (see Appendix [U] for numeri- 
cal estimates), the response of the spinon-gauge field sys- 
tem changes dramatically at low temperatures such that 
the Landau quantization of the spinons in the static inter- 
nal field is not smeared by the temperature. These effects 
are similar to strong electronic magnetism in quantizing 
field at low temperatures familiar in magneto-oscillation 
studies of metals^iiSii 3 . (which we also review in Ap- 
pendix El- in particular, the homogeneous state with 
continuously varying internal field becomes unstable be- 
low some temperature, which for the k-(ET) 2 Cu2(CN) 3 
spin liquid we estimate to be several Kelvin in typical 
laboratory fields. The instability regime is significantly 
wider in the spinon-gauge system than in metals because 
the internal gauge field stiffness is so much smaller than 
that of the physical electromagnetic field. Crudely, the 
spinon states with integer Landau level filling are more 
stable than the states with a continuously varying filling, 
and it becomes advantageous for the internal gauge field 
to adjust itself discontinuously to achieve this. This in- 
stability also preempts the possibility of a direct observa- 
tion of the spinon Fermi surface using magneto-oscillation 
probes, but can be viewed as an extreme manifestation 
of such magneto-oscillations. 

The objective of the present paper is to characterize 
the above two predictions in the spinon-gauge system fo- 
cusing on the k-(ET)2Cu2(CN)3 material. The paper is 
organized as follows. In Sec. UTI fand also Appendix 0, 
we describe the effective spin Hamiltonian derived from 
the triangular lattice Hubbard model in the magnetic 
field. In the main part of the paper, Sec. IIIII we per- 
form a spinon mean field study of the effective Hamil- 
tonian. We start with a review of the system in zero 
field^ and progressively add the needed ingredients for 
the discussion of the response to the magnetic field such 
as the appearance of the static internal gauge field and 
the importance of the spinon Landau level quantization 
at low temperatures. We then argue that a homogeneous 
spin liquid state with continuously varying static internal 
flux becomes unstable at low but experimentally relevant 
temperatures of few Kelvin. We finalize this section with 
a more careful discussion of the physical setting in the 
real system. Helpful connections with the gauge theory 
description and analogies with magnetic interaction ef- 
fects in metals are summarized in Appendixes IU1 and IDl 
Throughout, we use the k-(ET) 2 Cu2(CN) 3 parameters as 
a guide for the relevant questions. In Sec. IIV1 we collect 
the main estimates for the k-(ET)2Cu2(CN)3 material 
and consider some experimental aspects. 



II. RING EXCHANGE HAMILTONIAN IN THE 
PRESENCE OF MAGNETIC FIELD 

The approach adopted in Ref. 9 and pursued here is 
to focus on the spin degrees of freedom when describ- 
ing the insulating state. This is achieved by consider- 



ing an effective spin Hamiltonian obtained from the mi- 
croscopic Hubbard model by a canonical transformation 
that projects out the double occupancy. The importance 
of the charge fluctuations is retained in the form of more 
complicated multispin exchanges. The spin system is still 
more amenable to analysis because there is much less dis- 
proportion between the relevant energy scales and the 
couplings in the effective Hamiltonian, unlike in the orig- 
inal Hubbard model. 

The effective Hamiltonian to order t 4 /U 3 on the 
isotropic triangular lattice reads (see also Appendix |A"|) 

J 2 ^Pl2 + J 4 (^234 +H.C.) (1) 

+ j" J2 s * ■ s i + J "' E s * ■ s s ( 2 ) 

««» mm 

+ ® C A h E 1 ( P 123 - H.C.) . (3) 

A 

Here wc use multispin exchange operators defined as 

P\2...n ■ | £71, CT 2 , • ■ • , <J n ) — > |<7 n , <7l , . . . , £7„_l) . FortwO 

spins, this reduces to the familiar Heisenberg exchange, 
Pi2 = P12 = 2Si • S% + \. The first two lines give the ef- 
fective Hamiltonian in the absence of the magnetic field. 
As discussed in Ref. 0, important terms in this Hamil- 
tonian are the nearest-neighbor two-spin exchanges and 
the ring exchanges around the rhombi of the triangular 
lattice. The corresponding coupling constants are given 
in terms of the Hubbard model parameters as 
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The second- and third-neighbor Heisenberg exchanges 
J" = -16t 4 /C/ 3 and J'" = U 4 /U 3 , though nominally 
of the same order as the ring exchange coupling, were ar- 
gued to play a minor role because of the weak correlations 
between such further-neighbor spins in the magnetically 
disordered state. 

The last line in Eq. © shows a new term that appears 
at order t 3 /U 2 in the presence of the magnetic field. This 
term involves three-spin exchanges around the elemen- 
tary triangles and is proportional to the enclosed flux 
with the coupling constant 



_ 6i J 



(5) 



The dimensionless flux is 



eBA&/(hc), where B 
is the field and is the triangle area. The Hamilto- 
nian Eq. Q is written to linear order in the external 
flux, which is assumed to be small, <C 1. Each 

triangle is counted once and is traversed in the same di- 
rection. We can write the J3 term more explicitly as 
sm^ga) i (P123 — H.c), where ^ s the complex phase 
of the loop product <i2^23^3i of the electron hopping am- 
plitudes in the field. In general, each contribution from 
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an exchange path that encloses flux is affected by the 
magnetic field (cf. Appendix [SJ. However, the couplings 
J2, J4, J", and J'", are modified only at quadratic order 
in $ ext . The exhibited three-spin terms represent the 
full effect linear in <& GXt to order t 4 /U 3 . The effective 
Hamiltonian contains only terms P12, (-P1234 + H.c), and 
* (P123 — H.c), which was achieved with the help of the 
following identities valid in the spin-1/2 case: 



P 123 + H.c. = P u + P 23 + P31 - 1 , 



(6) 



P1234 - H.c. = - (P 123 + P234 + P341 + P412 - H.c.)(.7) 

The three-spin operator i (P123 — H.c.) has a simple 
physical meaning - it represents the spin chirality: 



i(P 123 -H.c.) = -4S 1 -S 2 



(8) 



Thus, the external field couples linearly to the spin chi- 
rality and therefore induces such chirality density in the 
system. 

In the above discussion, we have not mentioned the 
original Zeeman coupling of the electrons to the magnetic 
field. Since the Zeeman term involves the conserved total 
S z , we can simply reinstate it in the final effective spin 
Hamiltonian. 

For the k-(ET) 2 Cu 2 (CN) 3 , we use t = 55 meV and 
U/t — 8.2 to obtain J2 ~ 7 meV, J4/J2 — 0.3, and 
J3/J2 — 0.7; the Heisenberg exchange coupling is J = 
2J2 ~ 160 K, which is comparable with the estimate 
J ~ 250 K in Ref. [J 

To summarize, when charge fluctuations are signifi- 
cant, the effective spin Hamiltonian in the presence of 
the applied magnetic field is modified beyond the direct 
Zeeman term and contains a linear coupling of the spin 
chirality to the external field. This orbital effect is im- 
portant when considering the response of such insulator 
to the magnetic field. 



III. SPINON MEAN FIELD DESCRIPTION 
A. General setting 

In Ref. we considered the situation with no mag- 
netic field and argued that the four-spin ring exchanges 
stabilize the spin liquid state with spinon Fermi surface. 
This state can be viewed as a Gutzwiller projection of a 
fermionic mean field state in which spinons hop on the tri- 
angular lattice. We performed a direct variational wave 
function study, and also provided an intuitive mean field 
argument for the stabilization of this state. 

A similar mean field treatment is pursued here to un- 
derstand the effects of the magnetic field. To simplify the 
discussion, we will focus on the J2-J3-J4 terms in i? e ff- 
As mentioned earlier, the second- and third-neighbor 
Heisenberg exchanges are not expected to have signifi- 
cant effect. 



Each spin-1/2 is represented in terms of two "spinon" 
operators /•[• , / j. , with the occupancy constraints^ 

frcrfrcr = 1 ■ (9) 

In this representation, the J2-J3-J4 Hamiltonian reads 

links 

+ J*J2 [(/l t a/l/3)(/2 t /3 /27)(/3V3^)(/4V4 a )+H.C." 
rhombi 

+ J 3 M^Dl [(/ 1 t a/l/j)(/i/»r)(/i r /3a)-H.C 



triangles 



For clarity, we use the form of the </ 3 term valid for gen- 
eral $ oxt (cf . Appendix^) , even though we are interested 



in the case of small $g* ^ 1- 

We consider trial spinon hopping Hamiltoniansi£*i£ 

fft ri al = -E(C' n0 ^ a /r' a +H.C.)-^ Mr /t a / ra (10) 

(rr') r 

with general hopping amplitudes t^ 1 , 11011 = (t s ^" wn )* . The 
occupancy constraints are implemented on average by 
the appropriate chemical potentials fj, r . The mean field 
scheme that we use here is to evaluate the above exchange 
terms by contracting only the fermions with the same 
spin index. This approach can be justified in a large- 
N fermionic generalization, since any other contraction 
is down by a factor of 1/N. While the present N — 
2 is not large, the other contractions do not introduce 
qualitatively new terms, and this scheme is expected to 
capture the relevant physics rather well. The mean field 
energy is 



Pmf 



-4 ff2 J 2 ^ 1X12 I 2 (11) 

links 

I634 J A ^2 [X12X23X34X41 + C.C.] 



rhombi 



853^3 J2 Sm ( $< l23)(-^)[X12X23X31 -C.C. 
triangles 



Here \ r 



- //t 



(fj-'fr) are nn k expectation values evalu- 



ated for one spin species, and the powers-of-two numer- 
ical factors appear from the spin summations. We also 
introduce phenomenological renormalization factors 32, 
<74, and g 3 , to keep us aware of the schematic charac- 
ter of the mean field treatment. In a more quantitative 
treatment, these factors can be estimated by matching 
to numerical evaluations with the Gutzwiller-projected 
wave functions. This is done in Appendix [51 where we 
find that 172 — 1-7, 34 ~ 8, and <? 3 ~ 5, reproduce fairly 
well such direct trial wave function computations. 

The mean field energy is to be minimized over the 
trial spinon hopping amplitudes t^", non . We focus on 
the states that describe translationally invariant spin 
liquids. 15 These are the projected Fermi sea state with 
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real hopping amplitudes and the so-called flux states 
with complex t!? 1 ,* 10 " realizing nontrivial "internal" fluxes 
through the hopping loops. The flux can be either uni- 
form or have a staggered pattern. The projected Fermi 
sea state is a special case with zero flux. Here, we are 
primarily interested in the uniform flux states since these 
are natural candidates in the presence of the external 
magnetic field if one starts with the zero-flux state in the 
absence of the magnetic field. 

The mean field energy per site for a uniform flux state 

is 

e mf = -12g 2 J 2 |x0| 2 - 9634^4|X0| 4 cos(0i23 4 ) 

- 32 5 3 t /3|X0| 3 sin($- t 3 )sin(4 n 2 t 3) . (12) 

Here, Xrr' = \x<t>\ exp(za rr <) and </> lnt is the flux of the "in- 
ternal" gauge field a rr /. Note that the link expectation 
values Xrr' obtain the same flux pattern as the input am- 
plitudes t^) non . The gauge theory language is explained 
in Appendix [HI 

B. Review of the zero- field case 

Let us first consider $ cxt = following Ref. 0- We 
find that for small giJi/fafo) < 0.69, the lowest energy 
state has </>a = tt/2 flux through each triangle, while for 
larger ring exchanges the best state has zero flux (see 
Fig-O- The corresponding numerical values of \x<p\ for 
the half-filled triangular lattice are 

Ix*=*/al = 0.2002, |x*=o| = 0-1647 . (13) 

It is knownii that flux states have large absolute value 
\x<j>\ and therefore good Heisenberg energies. On the 
other hand, the ring exchanges are directly sensitive to 
the placket fluxes and dislike any fluxes, as can be seen 
from Eq. (|12(l . This is why the zero-flux state is stabi- 
lized for large J4 / J 2 ■ As explained in Ref. 0, we do not 
consider so-called dimer states even though formally for 
fixed 32 J2 these have the lowest mean field energy in the 
range 34^4/(32^2) < 2.4. Our reasoning is that J2 is 
fixed, and to go from the mean field evaluations to the 
projected wave functions, the renormalization factor 32 
is larger for the translationally invariant states. 

For later convenience, Fig. [21 shows the behavior of 
|X0 1 evaluated numerically for the triangular lattice flux 
states. It is useful to note the enveloping function for 
small fluxes, which we find to be 

|3Wl<IXo|(l+c&) (14) 

with c w 0.1. In particular, it follows that for small 
34^4/(32^2) < 0.56 - hatched region in Fig.Q~ the zero- 
flux state is also unstable against introducing weak fluxes 
(while the 7r/2-flux state is the global translationally- 
invariant minimum till a somewhat larger value of 0.69). 

A direct Gutzwiller wavefunction study gives that the 
zero-flux state has lower energy than the 7r/2-flux state 




7t/2-flux no flux (spinon FS) 
0.7 g4V(g2 J 2 ) 

FIG. 1: "Phase diagram" from the mean field energy op- 
timization over translationally invariant states in the ab- 
sence of the magnetic field. In the hatched region, the zero- 
flux state is also unstable against introducing small inter- 
nal flux. For the k-(ET)2Cu2(CN)3 parameters, we estimate 
g4.J4.Kg2 J2) * 1.2 - 1.5. 
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1/V = 2<VTC 

FIG. 2: Mean field numerical data for uniform flux states on 
the half-filled triangular lattice. The figure shows the abso- 
lute value of the bond expectation value |x<^| = \(fl/fr)\ as 
a function of the spinon filling factor v with respect to the 
internal flux (one spin species is considered). Fhe dotted line 
is the enveloping function Eq. 1141 , while the dash-dotted line 
is the full model function Eq. 1171 . The data plotted corre- 
spond to rather large </> lnt , but the same behavior is expected 
to continue for small fluxes. 

for J4 > 0.145J2. From this and similar comparisons we 
estimate that 34/32 is roughly 4 to 5 - see Appendix IbI 
for more details. 

In the k-(ET) 2 Cu2(CN) 3 compound, we have J4/J2 — 
0.3, which is not far from the regime where it would 
be advantageous to spontaneously generate such inter- 
nal flux. We therefore reason that the spinon Fermi sea 
state in the k-(ET)2Cu2(CN)3 compound is rather sus- 
ceptible to the internal flux generation. In particular, we 
expect some enhancement in the response to the external 
magnetic field, to which we now turn. 



C. Mean field over homogeneous flux states in the 
presence of the magnetic field 

Consider the mean field energy Eq. I|12|l with nonzero 
but small <j> GXt . Let us first try the enveloping function 
Eq. I|14|) for small 4> mt . Expanding e m f to quadratic order 



5 



in <j) mt and minimizing over the internal flux, we find 



53^3 



l2g 4 J 4 \xo\ 



1 - 2c 1 



92 J-2 

I6g<i -Ma (i 



4> c 



7 $ c 



(15) 

The different ingredients have the following physical ori- 
gin. The J3 in the numerator represents the coupling 
of the magnetic field to the internal gauge flux through 
the corresponding three-spin term. On the other hand, 
the J4 in the denominator represents the stiffness of the 
internal gauge field that originates from the four-spin ex- 
changes. For the k-(ET) 2 Cu2(CN) 3 parameters, we es- 
timate 33 ^3/(12(74 J4 1 Xo I) — 0.7. The term in the square 
brackets in the denominator reflects the enhanced suscep- 
tibility to the fluxes discussed above; it can be also viewed 
as a suppression of the internal gauge field stiffness. In 
particular, when this term goes to zero, the spinon Fermi 
sea state becomes unstable to spontaneous flux genera- 
tion even in the absence of the external field. Given the 
proximity of the competing states, it is therefore reason- 
able to estimate this number to be of order one-half, so 
for the k-(ET)2Cu2(CN)3 material we roughly estimate 



7k-(ET) 2 Cu 2 (CN) 3 



1 -2 . 



(16) 



Thus, we conclude that due to the triangular ring ex- 
changes and the proximity to the flux (or antiferromag- 
netic) instability, the effective orbital field seen by the 
spinons is comparable and can even be larger than the 
applied magnetic field! 

We now examine the mean field energy Eq. i|12|) more 
carefully. In the presence of the static internal gauge 
flux, the spinon spectrum consists of Landau bands. For 
either spin species, the filling factor of these Landau lev- 
els is v = 7i"/ (20^*), and the flux states are special when 
v is integer. Thus, the bond expectation value \x<t>\ has 
upward-pointing cusps at these fluxes as can be seen in 
Fig. [3 The states with integer v are therefore expected 
to be more stable. In fact, when the mean field energy 
is minimized with respect to <p lnt using the material pa- 
rameters quoted above, we find that the optimal such flux 
does not change continuously but instead goes through 
this discrete set corresponding to the integer filling of 
the spinon Landau levels. This is shown in Fig. The 
overall magnitude is still given by Eq. (fTBT) . 

Let us exhibit the structure of the discussed mean field 
in a more transparent way. The reason behind the ele- 
mentary manipulations below is to draw a connection 
with the so-called strong magnetic interaction effects in 
the magneto-oscillation studies in metals. These are re- 
viewed in Appendix [D] which contains useful interpreta- 
tions of the expressions below. At the expense of being 
repetitious, what follows will allow us to estimate appro- 
priate effective parameters, which we will then use in our 
discussion of the spinon system at finite temperature. 

First of all, we model the behavior of jx^l, Fig. [3 as 



\X<j>\ 



IXo|(l 



■c0a) 



" osc 



(17) 



20 



15 



-a- 

CM 



Mean field energy minimization 

g 2 J 2 =1.7 
g 3 J 3 = 3.5 

g 4 J 4 = 2.4 



10 



ii 

_ 

> 



(7=1-5) 



10 

V 



15 

ext 71/(2*?") 



20 



25 



30 



FIG. 3: Result of the mean field energy minimization with 
respect to <jf> mt in the presence of the magnetic field at zero 
temperature. The optimal spinon filling factor 1/ is plotted 
vs the nominal filling factor v cyit for the applied magnetic 
field. The figure shows a sequence of first-order transitions 
with i/ nt stepping through integer values. The parameters 
roughly correspond to our estimates for the K-(ET)2Cu2(CN)3 
from Sec. [H] with appropriate renormalization factors from 
Appendix[5] The dotted line shows the overall trend Eq. 115^ . 
The data plotted correspond to unrealistically large fields, 
but the same behavior is expected to continue for smaller 
laboratory fields. 



where we separate the envelope Eq. Q14[l. which is the 
band structure effect, and the oscillating part S^ sc due to 
the discreteness of the spinon Landau levels. The latter 
is well approximated by the familiar expression for the 
oscillating piece of the energy of the continuum Landau 
problem, 



[y - k){k + 1-1/) 



(18) 



for the filling factor v between integers k and k+1. This 
approximation is also plotted in Fig. [21 with a suitably 
chosen numerical amplitude for S^ 50 . In particular, from 
the numerical data we estimate the cyclotron mass for 
the half-filled triangular lattice band to be m c k, t^p; non 
in units with the lattice spacing set to one, while the 
corresponding cyclotron frequency at filling v is 



fku c w 3.57 t s 



(19) 



Returning to the mean-field energy Eq. 1)12(1 . we ex- 
pand to quadratic order in </> mt bearing in mind that the 
oscillating piece Sjl sc has an overall magnitude (cV nt ) 2 : 



emf = K4 [(^ t ) 2 -2 7 $ e A x VA] +■ 



«4 y 



— 1 9 + ^OSC 
1Z t S pj non 



with the effective parameters 

k 4 = 192g 4 J 4 \xQ\ 4 (l-2cy) , 

Spinon = 2g 2 J2|Xo| + 32g 4 J 4 |xo| 3 



(20) 
(21) 



(22) 
(23) 
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Here 7 is defined in Eq. (|TB|) : y = 1 +52^2/(1634 J4IX0I 2 ); 
and y' = 2y/ (1 — 2cy). We remind that the energy here 
is per triangular lattice site, and see that K4 is the cor- 
responding measure of the stiffness of the internal gauge 
field — see Appendix |0 We also see explicitly how the 
magnetic field acts to induce the internal gauge flux. On 
the other hand, using Eq. (|17|l , we can interpret e osc as an 
oscillating piece of some fermion kinetic energy of spinful 
fermions hopping on the triangular lattice with ampli- 
tude 4pinon- It is useful to remember that we are ana- 
lyzing the mean field energy Eq. H12[l . and from this we 
are separating out what looks like an effective spinon ki- 
netic energy. The division of the mean field energy into 
the gauge field and the spinon parts is somewhat arbi- 
trary, but represents a useful separation by the character 
of their dependence on <fi mt . In particular, we can now 
see the similarity with Eq. 1)03(1 and thus connect with 
the studies of magnetic interaction effects in metals sum- 
marized in Appendix iDl The above equations constitute 
the main result of this section. 

To emphasize the variational energetics character of 
the mean field procedure and the intrinsically common 
origin of the gauge field and the spinon parts, we do not 
use such suggestive separation explicitly in the treatment 
below. A formal condition for an extremum of e m f reads 



7*2* = 



y 



2|xo| 



6£ 



4a] 



(24) 



which can be solved graphically by plotting the right- 
hand side as a function of </> mt and seeking crossings 
with the horizontal line 7$ cxt . This can be compared 
with Eq. I|D4|) in Appendix IDl Here we remark that the 
oscillations of the "magnetization" 9Jt osc have roughly 
the sawtooth pattern familiar in two dimensions, with 
the amplitude which can be estimated from the data in 
Fig. |21 Note that the internal gauge field stiffness is im- 
plicitly included in 97t osc , cf. Eq. (|D4(1 . We estimate that 
the corresponding |9Jt osc | ma x is as large as 1. In partic- 
ular, we conclude that at zero temperature essentially 
each oscillation period supplies such graphical solutions 
of Eq. I|24l) . Solutions are found even in the region with 
v ~ 1 (corresponding to very large internal flux) even 
when $ ext is small. Remembering the cusps in this 
means that each integer filling gives a locally stable min- 
imum, which we can also verify directly by plotting e m f 
as a function of <p lnt . The global minimum is given by an 
integer filling near that corresponding to Eq. 1)150. This 
provides a more complete description of our direct mini- 
mization result Fig. |2| at zero temperature, and is useful 
when generalizing to finite temperature. 

The multiple solutions signify an instability of the flux 
states with continuously varying <fi lnt . The instability oc- 
curs whenever the derivative of the right-hand side of 
Eq. I|24|) becomes negative: 



dWl° 



> 1 



> 2 |2H osc | max > 1 



(25) 



This is just a more formal restatement of our discussion 
in the preceding paragraph: given |9Jt osc | max ~ 1, the 
instability condition is satisfied all the way to very large 
internal fields corresponding to v ~ 1. This is unlike the 
ordinary metals where the largest magnetic field that al- 
lows such an instability is typically much smaller. As 
explained in Appendix [DJ the origin of this is the very 
small internal gauge field stiffness, which is implicitly 
built into the numerical estimate |9Jt osc | max ~ 1 in the 
above formalism. However, such large numerical value of 
|9Jl osc | max is unavoidable since it is the same mean field 
energetics that stabilizes the Fermi sea state and deter- 
mines both the gauge field stiffness and the spinon kinetic 
energy. 



D. Mean field description at finite temperature 

The preceding analysis was performed at zero temper- 
ature. Also, we ignored the Zeeman spin coupling to 
the magnetic field. We first note that for integer fill- 
ing v the spinon spectrum has a gap given roughly by 
the corresponding cyclotron frequency. Such an integer 
quantum Hall state of spinons is in fact a chiral spin 
liquid 18 induced by the magnetic field, and the above 
treatment suggests a sequence of first-order transitions 
stepping through such states rather than a continuous 
variation of 4> lnt (see Fig. |3J). Finite temperature smears 
the effects of the discreteness of the Landau levels, and 
the instability of the flux states with continuously vary- 
ing 4> lnt becomes weaker. 

To proceed more quantitatively, we again use mean 
field and apply standard magneto-oscillation results at 
finite-temperature in two dimensions supplemented with 
our estimates of the effective parameters such as the cy- 
clotron frequency Eq. JT5J. Thus, we can roughly in- 
corporate the effects of the finite temperature and the 
Zeeman spin splitting by multiplying the oscillating mag- 
netization 97t osc by the corresponding suppression factors 
Rt and Rs from Eq. (|D2|I in Appendix IDl The instabil- 
ity condition Eq. (|25|l then reads 



R T R S v 2 \Tt° 



> 1 



(26) 



For a fixed external magnetic field, we can use Eq. i|15|) 
to estimate the corresponding spinon filling v in the in- 
ternal gauge field. Focusing on the temperature suppres- 
sion factor, the homogeneous flux state becomes unstable 
when the temperature is below 



r,, 



stab. 



xfnO c 
27T 2 



sinh(a;) 



i/|3Jt° 



(27) 



For high filling factors v and using 97t osc | max ~ 1, we 
estimate with logarithmic accuracy x ~ 2 log(^). We give 
numerical estimates for the k-(ET)2Cu2(CN)3 material 
in Sec. US 

Let us briefly think what happens beyond the mean 
field. The chiral spin liquid states are stable topologi- 
cally ordered phases^ at zero temperature, but strictly 
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speaking do not survive on long length scales at any 
nonzero temperature in two dimensions. Instead, the 
spin system can be continuously connected to the fea- 
tureless high-temperature paramagnet. In this case, we 
interpret the mean field estimate Eq. I|27[) as delineating 
the temperature-field regime in which the spin system 
becomes particularly "soft" in its response to the mag- 
netic field. It is also possible for a T = first-order 
transition between chiral states to extend to small finite 
temperatures even though the paramagnetic phases on 
both sides do not have a symmetry or topological dis- 
tinction anymore. Note also that in the above schematic 
treatment, we have not considered the temperature de- 
pendence of the effective parameters For example, in the 
formal mean field, the effective i s P inon decreases with 
temperature and disappears above T ~ 50 — I00K for 
the k-(ET)2Cu2(CN)3 parameters (this can be viewed as 
a crude estimate of the temperature below which quan- 
tum spin correlations develop). Such detailed questions 
on the behavior of the spin system at finite temperature 
remain open for future investigations. 



E. Inhomogeneous state in the physical system at 
low temperatures 

We conclude the mean field description with a dis- 
cussion of some experimentally relevant aspects. Our 
analysis indicates the instability of the homogeneous flux 
states with continuously varying (f> lnt at temperatures be- 
low Tinstab.- Throughout, we treat the external magnetic 
field as fixed, and the instability represents a strong non- 
linear back action of the spinons onto the internal gauge 
field in the presence of such fixed external source. Treat- 
ing the magnetic field as fixed is justified since, as we de- 
scribe in Appendixes IT?! and Hfl the electromagnetic field 
is six to seven orders of magnitude more stiff than the 
internal gauge field and will not adjust itself to the elec- 
tronic system until much lower temperatures. Of course, 
there is a slight effect on the local magnetic field depend- 
ing on the state of the spin system, but for example for 
the uniform flux mean field state we estimate the or- 
bital contribution to be one order of magnitude smaller 
than the Pauli spin contribution^ The physical setting 
therefore has fixed <I> cxt . In particular, we do not have 
access to different sample "demagnetization geometries" 
considered for the magnetic interaction phenomenon in 
conventional metalsAi*^ - the internal gauge field does 
not "leak out" of the sample. 

In such setting in an ideal crystal system at zero 
temperature, we predict a sequence of first-order 
transitions j2£ but each phase is still a uniform flux state 
with the corresponding integer v. This situation, how- 
ever, is highly susceptible to large-scale imperfections of a 
real system, and an outcome with many domains is likely. 
Another important consideration is the possible crystal 
mosaic in the sample, since the discussed orbital effects 
are determined by the component of the magnetic field 



H±_ normal to the two-dimensional plane. A detailed 
characterization of the possible inhomogeneities clearly 
requires much more material knowledge. 



IV. APPLICATION TO k-(ET) 2 CU 2 (CN) 3 

Throughout, we have used the k-(ET) 2 Cu2(CN) 3 pa- 
rameters to motivate various approximations and make 
the otherwise formal discussion more physical. Here we 
want to focus more on the material itself and collect the 
relevant numbers in one place. 

We first want to point out that the material is rather 
clean. Thus, we estimate UfI ^ 50 in the metallic phase 
at 0.8 GPa. Shubnikov-deHaas (SdH) oscillations are 
observed 21 at temperatures around 1.5 K, which is an- 
other indication of the material quality. The SdH signal 
is consistent with the Fermi surface and the one-band 
model description. For reference, the Landau level fill- 
ing factor for the k-(ET)2Cu2(CN)3 material is given 
by i/ cxt = 3595/#[Tesla]. In Appendix El we esti- 
mate the characteristic magnetic field Ho ks 5 T and 
temperature T^ m w 0.16 K for the magnetic interaction 
instability22i22i2£ in this metallic phase coupled to the 
electromagnetic field. 

Turning to the spin liquid phase, in Ref. we used 
the low-temperature susceptibility to estimate the spinon 
hopping amplitude i S pinon — 350 K. We can also use the 
mean field estimate Eq. I|23l) . which gives a smaller but 
reasonable value t sp i non ~ 100 K. Using Eqs. (|T§|| and 
(HZJ), and u {at = v cxt /j with 7 ~ 1 - 2, we obtain the 
characteristic temperature at B = 8 T to be Tinstab. — 
1 — 2 K. Above this temperature, we expect uniform flux 
stated 

We now summarize qualitative predictions for the in- 
homogeneous state induced by large magnetic fields at 
low temperatures below Ti nsta b- As discussed in the para- 
graph following Eq. 124fl . the characteristic field Hq for 
the spinon system is beyond any practical fields, and 
we predict that in the laboratory fields the instability 
temperature will increase roughly linearly with the ap- 
plied field, Eq. I|27|) . Manifestations of this instability 
become more pronounced with increasing field and de- 
creasing temperature. Since this is an orbital effect, we 
expect distinction for fields perpendicular and parallel to 
the conducting planes. One consequence of this scenario 
is that magneto-oscillation measurements in their usual 
sense cannot be used to detect the spinon Fermi surface 
in the spin liquid state considered here. This is because 
the regime where such oscillations can become visible is 
likely preempted by the instability of the homogeneous 
state. Since the instability is in some sense an extreme 
manifestation of the magneto-oscillations, conditions of 
the magnetic field uniformity and the crystal mosaic are 
important experimental considerations, as discussed in 
Sec. IIII El Other measurements in the spin liquid phase 
performed in high magnetic field and at low temperature 
may also require more careful interpretation. Finally, 
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some different direct probes of the Fermi surface such as 
magneto-acoustic resonance that are less sensitive to the 
field homogeneity may still be possible. 

V. CONCLUSIONS 

The main content of this work was summarized in 
the introduction, and possible application to the k- 
(ET)2Cu2(CN)3 experiments including some numerical 
estimates was discussed in the preceding section. Here 
we want to conclude with a more general point. While 
there has been a significant progress in the theoretical 
understanding of exotic quantum phases, the lack of the 
material evidence is viewed as a major obstacle. Re- 
cently, several candidate spin liquid systems have been 
found in frustrated spin systems ^S^S^ of which the k- 
(ET)2Cu2(CN)3 is a very promising example. This ma- 
terial motivated and guided the detailed theoretical con- 
siderations in the present work on the magnetic field re- 
sponse of the spin liquid state stabilized near the Mott 
transition. Without such a guide, the proposed effects 
would be hard to anticipate. More work matching ex- 
periment and theory focusing on material properties will 
likely be fruitful in developing our understanding and in 
further pursuits of such unusual quantum phases. 

One specific motivation for the present study has 
been the anomalous NMR line broadening reported in 
Refs. H and H While some of our expectations for the 
inhomogeneous state at low temperatures resemble the 
experimental phenomenology, the gradual development 
of the broadening starting from rather high temperatures 
and also the magnitude of the inhomogeneous fields are 
more suggestive of a stronger impurity mechanism. Ex- 
periments in parallel field can further clarify the role of 
the impurity and orbital effects. The main predictions 
in the present paper of the spinon orbital field and the 
fragility of the spinon Fermi sea state in the applied mag- 
netic field represent intrinsic response of the spin liquid 
and are expected to play an important role at low tem- 
peratures. It is hoped that our proposals will stimulate 
further experimental and theoretical questions. 

Note added: Recently, the author learned of Ref. [2^, 
which also discusses the coupling of the external field to 
the spin chirality and considers possible effects in some 



different contexts for Hubbard systems with triangles. 
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APPENDIX A: HUBBARD TO HEISENBERG IN 
THE PRESENCE OF MAGNETIC FIELD 

This appendix gives the effective spin Hamiltonian to 
order t 4 /U 3 for a general lattice Hubbard system with 
complex hopping t rr r — \t rr i \e lArr ' , which extends pre- 
viously available results^ to the case in the presence of 
the magnetic field. 

The starting point is the Hubbard Hamiltonian at half- 
filling 

H = U^2 »M»M - £ trr'cXaCr'a , (Al) 
r rr' 

with t r ' r = t* r> . The hopping part is treated as a per- 
turbation, and a canonical transformation is performed 
into the sector with no double occupancy. The effective 
Hamiltonian to order i 4 /£/ 3 reads 

H cS = H 2 + H 3 + H 4 , (A2) 

= £ ^(P 12 - 1) , (A3) 
•— • 

6Im(ti 2 t23*3i) . , D tj x i \ a\ 

H * = Jj2 1 ( P 123 - H - C I ( A4 ) 

6 



H 4 



•-• 1,2,3 

ERe(il 2 £23i34^4l) r on / D 
— IP — L 20(P] 



u 3 

1234 + H.C.) ' 



(P 2 3 - 1) 

- 12(P 12 + P 23 + P34 + P41 + P13 + P24) + 32 

I 



(A5) 
(A6) 



Here Pi 2 , -P123, and P1234 denote two-spin, three-spin, and four-spin exchange operators respectively. The latter 
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two operators move the spins in a ringlike manner. For 
example, the three-spin exchange in terms of the fermions 
is P123 = (cla c i/3)( c 2/3 c 27)( c 37 c 3a) ; and acts on the spins 
as P123 : o"2, C3) — > \<J3, ai, a%). To simplify the final 
expressions, we used the identities Eqs. (|6l7fl specific to 
the spin-1/2 case. 

The sum in H2 is over all bonds of the lattice. The 
sum in H3 is over all three-site loops, where each group of 
three sites connected by the links enters the sum precisely 
once. Similarly the sum in the second line for H4 is over 
four-site loops, which are counted with no base point 
or orientation. Finally, the primed sum in the first line 
for H4 is over distinct sites 1,2,3 with nonzero £12 and 
ti3 (each contribution will therefore appear two times). 
Eq. © in the main text is obtained by specializing to the 
isotropic triangular lattice. 

Note that the complex phase of a loop product of t rr > 
measures the flux of the magnetic field through the loop. 
Without such fluxes, H3 vanishes and after some trans- 
formations we reproduce the result of Ref. |27l 

As a side remark, we observe that in the absence of 
three-site loops, there is no linear coupling to the mag- 
netic field to this order in t/U. One can also show that for 
a half-filled Hubbard model on a bipartite lattice, the ef- 
fective Hamiltonian in the presence of the magnetic field 
can only contain terms that are even in the applied field. 

Finally, it is also useful to give an expression for an 
effective local electrical current operator. In the original 
Hubbard model, the current on a link 1 — ^ 2 is Ii^2 = 
■§l(ti2c\ a C2a — H.c). Upon the canonical transformation 
to the lowest nontrivial order in t/U , the current operator 
becomes in terms of the spin variables 



rcff e 6Re(ii2^23i3i) . 

^^iX TP 1 

3 



P123 - H.c.) . (A7) 



We can use this expression e.g. to estimate the physical 
magnetic field produced by the underlying orbital motion 
in the chiral spin states discussed in the main text^S 



APPENDIX B: RENORMALIZED MEAN FIELD 
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FIG. 4: Expectation value of the two-spin exchange P12 on 
each link evaluated in the Gutzwiller-projected uniform flux 
state. The data is for commensurate fluxes on triangular lat- 
tice cuts with sizes 36 x 12 (squares) and 48 x 12 (circles) with 
periodic boundary conditions. The line shows the renormal- 
ized mean field estimate with = 1.7 (the mean field data 
is obtained using Fig.|2j. Since we are interested only in the 
variation with c/> lnt , a constant offset is added to the mean 
field values. 
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1/V = 2(^/71 

FIG. 5: Same as in Fig. 2] but for the ring exchange operator 
A234 + H.c. around a rhombus. 



In this appendix, we estimate the renormalization fac- 
tors (72 j 54, and <?3 used in the mean field treatment in 
Sec. IIIII We consider the uniform flux states and mea- 
sure the relevant expectation values of P12, (P1234+H.C.), 
and i(P\2z — H.c.) in the Gutzwiller-projected wavefunc- 
tions. The variation of the expectation values with c/> mt is 
matched with the corresponding mean field estimates — 
see Eq. (|11|) . Figures 01 [SI and show this matching and 
are self-explanatory. We note that the Gutzwillcr wave- 
function evaluations can be performed only for small sys- 
tem sizes and large values of </> mt . However, we expect 
the same trend to remain for small fluxes as well. For 
small fluxes, an analytical mean field treatment is then 
performed (Sec. lIIl|l using the estimated renormalization 



factors. This is the main idea of the so-called renormal- 
ized mean field approach^ 



APPENDIX C: SPINON-GAUGE THEORY 

In the main text, we often use the gauge theory 
language^ and it is useful to state the connection more 
explicitly. We follow the treatment of Ref. Kj as an 
example. One starts with a slave-rotor representation 
c ra = fra el6r °f the Hubbard model and obtains the fol- 
lowing effective theory for the spinons f ra and chargons 
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FIG. 6: Same as in Fig. [I] but for the three-spin operator 
i(Pi23 — H.c.) around an elementary triangle. The orientation 
of the triangular loop 1 — > 2 — > 3 —* 1 coincides with that of 

/int jint 
<Pl23 = <PA ■ 



<i(P 123 -h.c.)> 



J dr{C f +Ce) 



(CI) 



5 



i(e r -e,.,-a,. r ,+A ox 5) 



The imaginary time is continuous, while the space lat- 
tice is retained to indicate similarity of the spinon part 
with that in Eq. H10[l. The spinons and the chargons are 
coupled with opposite charge to the internal gauge field 
(a, a ), whose spatial components represent fluctuations 
of the phase of the spinon hopping field, t s 7 p r ) non — tfe larr ' , 
while the temporal a® implement the slave-particle con- 
straints. The coupling to the external field A cxt is in- 
cluded schematically by assigning the electron charge to 
the chargon field. 

Assuming the chargons are gapped and integrating 
them out, we are left with the spinon-gauge system where 
the gauge field now has some stiffness K, 



cxt\2 



£„ = -(VAa-VAA m ) 



0\2 



(C2) 



Here VAa = d x a y —d y a x = b; we separated this "spatial" 
part for later convenience. 

Our mean field treatment Eq. (|ll(l roughly corresponds 
to considering the spinon-gauge system with static but 
possibly spatially varying gauge field a. Observe that 
the internal gauge field obtains its stiffness after integrat- 
ing out the massive chargons, and also observe how the 
magnetic field couples to the internal gauge field. This 
is similar to our discussion of the effective spin Hamil- 
tonian, where the crucial ingredient is also charge fluc- 
tuations. Note however that the effective Hamiltonian 
treatment is more explicit. Thus, the gauge field stiffness 



originates primarily from the four-spin ring exchanges, 
while the coupling to the magnetic field comes from the 
three-spin processes - see Eqs. (Illll and l)2()[l. In the 
schematic gauge theory derivation this distinction is not 
being made. The result j = 1 of such schematic treat- 
ment roughly agrees with our estimate using Eq. Ijl5|l . 
but this is rather fortuitous - e.g., as we mention to- 
wards the end of Appendix 1X1 in a Hubbard model on a 
lattice with no elementary triangles the linear coupling of 
the internal gauge flux to the applied magnetic field does 
not appear even to order t 4 /U 3 . To capture the detail 
present in our Eq. (|15l) in the gauge theory approach, one 
needs to go back to the original slave particle rewriting 
and carefully consider saddle point conditions and inte- 
grations over massive fields when deriving Eq. IjClf) in the 
presence of the external field. As an example, the sup- 
pression of the gauge field stiffness due to the proximity 
to the flux phase discussed after Eq. (|15fl originates from 
some modes which become soft near the flux phase. The 
effective spin Hamiltonian approach is more transparent 
in this respect and also allows quantitative estimates of 
the physical quantities. 

In particular, from Eq. 12211 we quote the internal 
gauge field stiffness in terms of the energy cost per tri- 
angular lattice site, e; n t = k^O^a*) 2 , 



k 4 ~ 192g 4 J 4 |xo| 



0.14 g 4 J 4 ~ 2 meV 



(C3) 



Here we used the Sec. UTI estimate J 4 ~ 2 meV, and used 
the Fig. estimate g 4 ~ 8. 

This can be compared with the corresponding stiffness 
of the physical electromagnetic (EM) field in the bulk of 
the k-(ET) 2 Cu 2 (CN) 3 compound, e EM = VB 2 /(8tt) = 



KEM«) 2 , 



8nA 2 



1.1 x 10 4 eV 



(C4) 



Here V = 850.6 A 3 is the 3d volume per triangular lattice 
site, j4a = 28.76 A 2 is the area of an elementary triangle, 
and = eBA/\/ (he) is the appropriate dimensionless 
external flux. 

We therefore conclude that the bare internal gauge 
field stiffness is ~ 10 7 times smaller than the electro- 
magnetic field stiffness. As we discuss in the following 
Appendix IdI this makes a dramatic difference in the re- 
sponse of the spinon-gauge system to the external field 
compared with that of electrons in a conventional metal. 



APPENDIX D: ANALOGY WITH MAGNETIC 
INTERACTION EFFECTS IN METALS 

The purpose of this appendix is to clarify the status of 
our mean field treatment in Sec. lIIII bv pointing out the 
analogy with the so-called magnetic interaction effects in 
the study of magneto-oscillations in metals. This mate- 
rial is available in textbooks. 11,12 The presentation below 
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parallels some of our discussion in Sec. Illll in a different 
language, but we hope that the reader will benefit from 
this duplication. 

For ease of reference, we first write down the oscillating 
part of the magnetization of a two-dimensional electron 
gas in a field H at a finite temperature T: 

M™(H) = R T R s ^sm(^pl . (Dl) 
irm c \ H J 

Only the main harmonic is shown. The magnetization 
is given per unit area, and n is the 2d electron density 
(including spin); m* is the cyclotron mass; F = nnch/e 
so that F/H — v is the Landau level filling for each spin 
species. The suppression factors Rt and Rs due to finite 
temperature and Zeeman spin splitting are 



Rt — 



2TT 2 T/hlO c 

sinh(27r 2 T/ hio c 



grnn 
2m e 



(D2) 



where g is the spin g-factor. 

The essence of the magnetic interaction effects is the 
back action of the electrons onto the electromagnetic 
field, which can be significantly enhanced by the oscillat- 
ing character of Ai. Specifically, consider fermions cou- 
pled to a dynamical gauge field a in the presence of the 
external field A ext , with the Lagrangian density Cf + C a 
given by Eqs. (|C1IC2(1 . The microscopic magnetic field 
is b = V x a and the external field is H ext = V x A cxt . 
In the case of electrons, a represents the physical elec- 
tromagnetic field and k is the appropriate EM stiffness; 
e.g., in a 2d layered material the magnetic field energy 
per unit area is d*b 2 /8ir where d* is the spacing between 
the layers. 

The key observation is that electrons see the average 
microscopic field B because the typical Larmor orbits 
are large. The mean field treatment is then to assume a 
static but possibly spatially varying such field B[r) and 
solve the Landau problem for electrons moving in this 
field. The mean field functional including the energy of 
the gauge field i q 11 ! 12 ! 13 



fimf — ^ferm(-B) + ~ 0® — iJ cxt ) 2 



(D3) 



where f2f crm (-B) is the free energy density for fermions in 
the static field B{r). The mean field functional is to be 
minimized with respect to B{r). If we look for a uniform 
solution, we obtain 



H cxt = B + I^^H = B — —M(B) , 

K OB K 



(D4) 



where in the last equation we defined the magnetization 
density M(B). For a 2d electron gas, the oscillating piece 
of the magnetization is given in Eq. (|D1(I . and the mean 
field treatment is to replace H — » B. 

Eq. (|D4J| is to be solved for B. It is customary to 
plot the right hand side as a function of B and seek the 
intersection with the horizontal line _ff cxt . When this is 



unique, the corresponding uniform field is the sought for 
stable solution. However, M.(B) contains an oscillating 
piece, and if we find that 



1 DM 
k~8B 



> 1 



(D5) 



then the solution is no longer unique, which signals an 
instability. Note that while the amplitude of the oscilla- 
tions |A^| max (T = 0) = neh/(2m*c) is small in metals 
compared to typical kH, the amplitude of the "suscep- 
tibility" dM./dB obtains an additional factor ~ F/H 2 , 
and the instability condition can be satisfied at low tem- 
peratures and not too large fields. After simple transfor- 
mations, the instability condition reads 



R T Rs n |^ 
2m* 



> -H 2 

~ 2 



(D6) 



where kp is the Fermi wavevector. When this condition 
is satisfied, the gain in the fermion energy when the field 
B is adjusted to obtain integer Landau level filling over- 
weights the magnetic field energy cost. Since Rt and Rs 
do not exceed 1, the above equation sets the maximal 
value Hq for the instability to occur at zero tempera- 
ture. On the other hand, at a finite temperature, we also 
require the Landau levels to be resolved, which is deter- 
mined by Rt ■ Roughly, for fields of order Hq , we require 
2ir 2 T < huj c with uj c set by Hq. 

The 2d Landau problem for free electrons can be solved 
beyond the above single-harmonic treatment, and Ref.l22l 
contains such expressions for the phase boundary of the 
instability region. For example, the field Hq is deter- 
mined from 



h 2 k 2 F 

2m* 



kH 2 



8tt 



(D7) 



Ref. |22j also contains formulae directly appropriate for 
the layered organic materials. As a numerical applica- 
tion, consider the metallic phase of the k-(ET)2Cu2(CN)3 
compound obtained under pressure of 0.76 GPa, in which 
magneto-oscillations were reported in Ref. |2l]. Using 
tp ~ lOOmeV, the cyclotron mass ratio r\ c = m* /m e ~ 4, 
the inter-layer spacing d* = 14.84 A, and ignoring sam- 
ple demagnetization effects, we estimate Hq ss 5 T; the 
peak temperature on the phase boundary of the instabil- 
ity regime is estimated to be Tdm ~ 0.16K. 

Various aspects of what happens once the ho- 
mogeneous state becomes unstable are discussed in 
textbooks>ii*i£ There is also a growing recent literature 
on the observation of Condon domains 23 in metals, sec 
Ref. and references therein. 

The point that we want to make is that our mean field 
treatment in Sec. IIIII has the same character as the de- 
scribed treatment of the magnetic interaction effects in 
metals. Thus, our mean field energy Eq. (|20|) can be read- 
ily identified with Eq. (|D3|I . while the self-consistency 
condition Eq. (|24|l corresponds to Eq. (|D4|I . Clearly, the 
treatment of the magnetic interaction effects requires tak- 
ing into account the energy of the gauge field. This is 
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implicit in our mean field treatment of Sec. IIHI At this 
mean field stage, we do not need to disentangle the gauge 
field energy from the total energy and worry about the 
emergent nature of the gauge field, and this makes the 
procedure more simple. 

An important difference between the spinon-gauge sys- 
tem and the electrons in metals coupled to the EM field 
is the very small gauge field stiffness in the spinon case, 
as we estimated in Appendix El Because of this, the in- 
stability condition is satisfied more readily in the spinon- 
gauge system. Thus, because the internal gauge field 
stiffness is 10 6 — 10 7 times smaller while the kinetic pa- 
rameters of spinons are not far from those of electrons, 
the above estimate of Hq is to be multiplied by a factor 
of order 10 3 . This is an impractically large field, and the 



spinon-gauge system in the laboratory is in the regime 
of much smaller fields. In this regime, the temperature 
at which the instability occurs in a given field H can be 
estimated as 

_xhuj c sinh(a;) nh 2 k 2 F /m* 

T instab . - , — ^— « kR2 . (D8) 

This temperature depends on the stiffness only logarith- 
mically, but as numerical estimates in Sec. lIIII show. even 
under the logarithm the seven orders of magnitude dif- 
ference in the relevant stiffnesses produces a much wider 
temperature range for the instability in the spinon-gauge 
system. 
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